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Abstract 

A  recently  developed  K-x  model  for  near  wall  turbulent  flows  is  applied  to  two  severe  test  cas¬ 
es.  The  turbulent  flows  considered  include  the  incompressible  flat  plate  boundary  layer  with  ad¬ 
verse  pressure  gradients  and  incompressible  flow  past  a  backward  facing  step.  Calculations  are 
performed  for  this  two-equation  model  using  an  anisotropic  as  well  as  isotropic  eddy-viscosity. 

The  model  predictions  are  shown  to  compare  quite  favorably  with  experimental  data 
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Introduction 


Two-equation  turbulence  models  have  become  increasingly  popular  for  the  calculation  of  prac¬ 
tical  aerodynamic  flows  which  can  have  important  applications  in  the  design  of  advanced  aircraft. 
The  reason  for  this  popularity  is  clear:  two-equation  models  incorporate  substantially  more  turbu¬ 
lence  physics  —  and  require  less  ad  hoc  empiricisms  —  than  the  older  algebraic  eddy-viscosity 
models  without  most  of  the  added  difficulties  that  arise  from  the  computational  implementation  of 
second-order  closure  models.  Among  the  existing  two-equation  models,  the  K-e  model  is  probably 
the  most  popular;  its  successes  and  limitations  have  been  fairly  well  documented  [1,2].  One  major 
difficulty  with  the  K-e  model  that  has  yet  to  be  fully  resolved  involves  its  application  to  near  wall 
turbulent  flows  where  inaccuracies  and  numerical  stiffness  problems  can  arise  (see  Patel,  Rodi  and 
Scheuerer  [3]  for  an  interesting  discussion  of  these  problems).  Recently,  Speziale,  Abid  and 
Anderson  [4]  developed  a  new  K-z  model  for  near  wall  turbulent  flows  wherein  modeled  transport 
equations  are  solved  for  the  turbulent  kinetic  energy  K  and  the  turbulent  time  scale  x.  While  the 
use  of  the  turbulent  time  scale  had  been  proposed  earlier  within  the  context  of  two  equation  models 
(cf.,  Reynolds  [5]  and  Zeierman  and  Wolfshtein  [6]),  none  of  these  previous  studies  rendered  a 
model  that  could  be  integrated  directly  to  a  solid  boundary  with  the  no-slip  condition  applied.  Spe¬ 
ziale  et  al.  [4]  demonstrated  that  their  new  K-z  model  yielded  excellent  results  —  and  was  compu¬ 
tationally  robust  —  when  integrated  directly  to  the  wall  in  an  incompressible  flat  plate  boundary 
layer  with  zero  pressure  gradient.  The  K-z  model  alleviates  the  problem  of  the  lack  of  natural 
boundary  conditions  for  the  dissipation  rate  in  the  K-e  model  since  the  turbulent  time  scale  x  van¬ 
ishes  identically  at  a  solid  boundary. 

The  purpose  of  the  present  paper  is  to  provide  a  more  comprehensive  testing  and  evaluation  of 
the  K-z  model.  Two  test  cases  are  chosen  to  evaluate  the  model:  (i)  the  incompressible  flat  plate 
boundary  layer  with  adverse  pressure  gradients,  and  (ii)  incompressible  flow  past  a  backward  fac¬ 
ing  step.  Unlike  in  the  previous  study  [4],  the  K-z  model  will  be  solved  with  an  anisotropic  as  well 
as  isotropic  eddy  viscosity.  It  has  been  demonstrated  during  the  past  decade  that,  in  some  cases, 
the  predictions  of  two-equation  models  can  be  enhanced  considerably  by  the  use  of  an  anisotropic 
eddy  viscosity  [7,8].  The  results  obtained  from  the  K-z  model  will  be  documented  in  detail  and 
recommendations  will  be  made  for  future  research. 


The  K-x  Model 

In  the  standard  form  of  the  K-z  model,  the  Reynolds  stress  tensor  x^.  =  uj  uj  (where  u\  is  the 
fluctuating  velocity)  is  of  the  form 
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=  3*Vvr 
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where  the  eddy  viscosity  is  given  by 


vr  =  C^Kx 
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and  u,-  is  the  mean  velocity,  K  a  ^  u  J  u-  is  the  turbulent  kinetic  energy,  e  =  v  du-  /dxj  du I  /dxj 

is  the  turbulent  dissipation  rate,  x  =  AT/e  is  the  turbulent  time  scale,  and  is  a  dimensionless  con¬ 
stant  that  assumes  a  value  of  0.09.  In  (2),/^  is  a  wall  damping  function  that  is  needed  for  asymp¬ 
totic  consistency  [4];  it  is  taken  to  be  of  the  form 

(3) 


where  Ret  =  —  is  the  turbulence  Reynolds  number  and  y+=  yu^N  is  the  usual  wall  coordinate 
given  that  itj  is  shear  velocity  and  v  is  the  kinematic  viscosity.  In  an  attached  boundary  layer,  the 
shear  velocity  is  usually  taken  to  be  Jx^ 7p  (where  t0  is  the  wall  shear  stress);  for  more  complicat- 

ed  turbulent  flows  with  separation,  the  shear  velocity  is  assumed  to  be  proportional  to  K  .  The 
wall  damping  function/^  goes  to  unity  sufficiently  far  from  the  wall.  In  the  anisotropic  version  of 
the  K-x  model,  nonlinear  corrections  are  added  to  the  eddy  viscosity  model  as  developed  by  Spe- 
ziale  [8].  This  leads  to  the  following  Reynolds  stress  model 

TV  =  \Kha - 2CA™u - icDcyiK^ (4V  + sv  -!%«.. )  (4) 

where 
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are,  respectively,  the  mean  rate  of  strain  tensor  and  its  frame-indifferent  Oldroyd  derivative  while 
Cq  is  a  dimensionless  constant  that  assumes  a  value  of  1 .68.  Of  course,  the  standard  eddy  viscosity 
model  (1)  is  recovered  in  the  limit  as  Cq  -*  0. 

The  Reynolds  stress  models  (1)  and  (4)  are  solved  in  conjunction  with  modeled  transport  equa¬ 
tions  for  K  and  x  which  take  the  form  [4]: 
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where 
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is  a  wall  damping  function  that  is  needed  for  asymptotic  consistency.  The  constants  Cej  and  Ce 2 
assume  the  values  of  1.44  and  1.83,  respectively,  whereas  for  their  preliminary  calculations,  Spe- 
ziale  et  al.  [4]  chose 


°tl  =  °x2  =  CK=  136  (10) 

with  the  understanding  that  these  constants  could  be  fine  tuned  to  a  non-equivalent  set  of  values 
when  a  future  optimization  over  a  variety  of  turbulent  flows  is  conducted.  When  aTi,  and  G% 
are  identically  equal,  this  K-x  model  reduces  to  an  equivalent  K- e  model  which  can  be  simply  ob¬ 
tained  by  replacing  x  with  K/e  in  (7)  and  replacing  the  x-transport  equation  (8)  with  the  equivalent 
equation 
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where  G£  is  equal  to  the  common  value  of  oT|,  cr^,  and  0%  given  in  (10). 
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Discussion  of  Results 

The  first  problem  that  we  will  consider  is  the  incompressible  flat  plate  boundary  layer  with  ad¬ 
verse  pressure  gradients  —  the  Samuel  and  Joubert  [9]  test  case.  This  represents  a  rather  severe 
test  since  the  adverse  pressure  gradient  is  strong  enough  so  that  the  flow  is  not  that  far  removed 
from  separation.  In  Figure  1(a),  the  mean  velocity  protiles  predicted  by  the  standard  K-x  model 
using  the  eddy-viscosity  representation  (1)  are  compared  with  the  experimental  data  [9]  at  two  sta¬ 
tions:  x-  1.76  m  and  x  =  2.87  m  (this  corresponds  to  free-stream  Reynolds  numbers  Rex  in  the 

range  of  2X106  to  4.X106).  It  is  clear  from  these  figures  that  the  model  compares  quite  favorably 
with  the  experimental  results.  In  Figure  1(b),  the  skin  friction  predicted  by  the  K-x  model  is  com¬ 
pared  with  the  experimental  data.  Again,  the  agreement  between  the  model  predictions  and  the  ex¬ 
perimental  data  is  extremely  good.  Unlike  most  other  two-equation  models,  the  K-x  model  does 
not  overpredict  the  skin  friction  for  this  adverse  pressure  gradient  test  case  (see,  Rodi  and  Scheuer- 
er  [10]  for  an  interesting  discussion  of  this  issue).  It  should  be  noted  that  the  use  of  an  anisotropic 
eddy-viscosity  does  not  yield  a  significant  change  in  the  results  for  this  case. 
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The  second  problem  to  be  considered  is  incompressible  turbulent  flow  past  a  backward-facing 
step  —  the  same  test  case  considered  by  Driver  and  Seegmiller  [1 1]  in  their  recent  experiments. 
For  this  flow  configuration  the  aspect  ratio  HjlH  is  8  and  the  Reynolds  number  based  on  the  inlet 
free  stream  velocity  and  step  height  is  approximately  38,500  (see  Figure  2).  The  calculations  were 
conducted  on  a  200x100  nonuniform  mesh  using  a  finite  volume  method  (see  Thangam  and  Hur 
[12]).  Here,  we  use  the  anisotropic  eddy- viscosity  model  (4)  since  it  has  been  recently  demonstrat¬ 
ed  that  such  models  give  rise  to  improved  predictions  in  recirculating  flows  [8,12].  In  Figure  3(a), 
the  computed  streamlines  obtained  from  the  anisotropic  K-x  model  are  shown.  They  indicate  re- 
attachment  at  approximately  6  step  heights  downstream  of  separation  —  a  result  that  is  in  agree¬ 
ment  with  the  experimental  data  [11].  In  Figure  3(b),  the  mean  velocity  profiles  obtained  from  the 
anisotropic  K-x  model  are  compared  with  the  experimental  results;  again  the  agreement  is  extreme¬ 
ly  good.  In  Figure  4(a),  the  profiles  of  the  turbulence  intensity  obtained  from  the  anisotropic  K-x 
model  are  compared  with  the  experimental  data.  While  there  are  some  discrepancies  between  these 
results  in  the  separation  and  recovery  zone,  on  balance,  the  agreement  is  fairly  good.  In  Figure 
4(b),  the  profiles  of  the  Reynolds  shear  stress  obtained  from  the  anisotropic  K-x  model  are  com¬ 
pared  with  the  experimental  data  [11].  Again,  while  there  are  some  discrepancies  in  the  results 
within  the  separation  and  recovery  zone,  the  model  predictions  are  still  satisfactory.  There  is  no 
question  that  this  near  wall  model  does  a  substantially  better  job  than  the  standard  K-e  model  which 
underpredicts  the  reattachment  point  by  an  amount  of  the  order  of  25-30%. 


Conclusion 

The  near  wall  K-x  model  of  Speziale  et  al.  [4]  has  been  tested  in  two  applications  involving 
incompressible  turbulent  flows  and  the  results  obtained  are  quite  encouraging.  The  standard  K-x 
model  —  with  an  isotropic  eddy  viscosity  —  yielded  excellent  results  for  the  challenging  incom¬ 
pressible  adverse  pressure  gradient  case  of  Samuel-Joubert  [9].  For  the  separated  flr  v  considered 
—  namely,  turbulent  flow  past  a  backward  facing  step  —  excellent  results  were  obtained  from  the 
K-x  model  based  on  an  anisotropic  eddy-viscosity.  In  fact,  we  believe  that  these  may  be  the  most 
accurate  calculations  for  this  backstep  problem  that  have  yet  been  obtained  from  a  two-equation 
model  without  the  ad  hoc  adjustment  of  constants. 

Future  calculations  are  planned  for  supersonic  boundary  layers  with  streamline  curvature  and 
separation  (shock-induced  or  otherwise).  For  these  flows,  the  K-x  model  will  most  probably  have 
to  be  applied  with  the  anisotropic  eddy  viscosity  model  and  it  may  become  necessary  to  incorporate 
dilatational  effects  into  the  model.  After  these  calculations  are  completed,  and  refinements  are 
made  accordingly,  we  will  have  a  much  better  idea  of  the  range  of  applicability  of  the  K-x  model. 
Nonetheless,  at  this  point  the  prospects  appear  to  be  quite  promising. 
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Figure  1.  Comparison  of  the  results  of  the  K<  Model  and  experiments  for  the  flat 
plate  boundary  layer  with  an  adverse  pressure  gradient  ( —  K<  Model, 
o  experiments  of  Samuel  &  Jonbert  [9]) : 

(a)  mean  velocity,  and  (b)  skin  friction  coefficient 
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Figure  2.  Schematic  of  the  diagram  for  turbulent  flow  past  a  backward  facing  step 
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Figure  3.  Results  obtained  from  the  anisotropic  /C-r  Model  for  incompressible 
turbulent  flow  past  a  backward  facing  step  with  ftyH  =  8,  Re  -  38,500.  (  K< 

Model,  o  experiments  of  Driver  &  Seegmiller  [11])  t 

(a)  mean  velocity  streamlines,  and  (b)  mean  velocity  profiles. 
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Figure  4.  Results  obtained  from  the  anisotropic  K-x  Model  for  incompressible 
turbulent  flow  past  a  backward  facing  step  with  ftyH  =  8,  Re  -  38,500.  ( —  /C-r 
Model,  o  experiments  of  Driver  &  Seegmiller  [11]): 

(a)  turbulence  intensity  profiles,  and  (b)  Reynolds  stress  profiles. 
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